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Abstract. We present observations of methanol lines in a sample of Class low mass protostars. Using a 1-D radiative transfer 
model, we derive the abundances in the envelopes. In two sources of the sample, the observations can only be reproduced by 
the model if the methanol abundance is enhanced by about two order of magnitude in the inner hot region of the envelope. Two 
other sources show similar jumps, although at a lower confidence level. The observations for the other three sources are well 
reproduced with a constant abundance, but the presence of a jump cannot be ruled out. The observed methanol abundances in 
the warm gas around low mass protostars are orders of magnitude higher than gas phase chemistry models predict. Hence, in 
agreement with other evidence, this suggests that the high methanol abundance reflects recent evaporation of ices due to the 
heating by the newly formed star. The observed abundance ratios of CH3OH, H2CO and CO are in good agreement with grain 
surface chemistry models. However, the absolute abundances are more difficult to reproduce and may indicate the presence of 
multiple ice components in these regions. 
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1. Introduction 

During the formation of a star, the gas undergoes impor- 
tant physical and chemical changes. In the prestellar phase, 
the gas is heavily depleted by accretion on grain mantles. 
When the gravitational collapse starts, the protostar warms 
the gas while ice mantle molecules are released into the gas 
phase. These released molecules trigger the formation of more 
complex m olecules through rapid reactions in the warm gas 
JCharnlev etalJl992h . 

While the importance of ice mantle evaporation and hot 
core chemistry around high mass protostars was known well 
over a decade ago, the existence of such regio ns in low mass 
protostar s has been established more recently. ICeccarelli et alJ 
<200QaE) have shown that H2O and H2CO abundances in the 
low mass protostar IRAS 16293-2422 are increased in the warm 
and dense inner part of the circumstellar envelope. In this re- 
gion, H2O and H2CO are evaporated from the grain mantles 
and are injected i n the gaseous phase. These findings were 
later confirmed by Sch oier et alJ J2002I) . who found strong ev- 
idence for the increase of several other molecular abundances 
(e.g. CH3OH, SO and SO2). Evaporation of H2O has also been 
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shown around one other Cla ss protostar, NG C1333-IRAS4 
dMaret et al.l2002l) . Recently. iMaret et alJ d2004l Paper I) car- 
ried out a survey of the emission of H2CO in a sample of nine 
low mass protostars, and concluded that in all the observed pro- 
tostars but one, H2CO is between two and three orders of mag- 
nitude more abundant in the inner part of the envelope than in 
the outer part, suggesting that the evaporation of grain man- 
tles is a common phenomenon in the inner parts of low mass 
protostars. 

While the importance of grain mantle evaporation is now 
well established in low mass protostars, the question of whether 

these molecules will form "daughter" molecules by hot core 

1 11 1 J 

chemistry (Charnl ev et alJU9 92) is still open. From a theoreti- 
cal point of view, if the protostar is undergoing collapse, only 
very rapidly formed second-generation c omple x molecules can 
be produced. As noted bv lSchoier et al 1 d2002h if the protostar 
is in free fall, the evaporated molecules will fall on the cen- 
tral object in a few hundreds years. Hence, while evaporation 
can still occur on such a rapid timescale, the composition of 
the gas phase is then 'frozen' in, because the dynamical time 
scale is much shorter than the time needed to fo rm second 
generation molecules (Rod gers & Charnlev! 120031) . However, 
observationally, the recent detections of complex O and N 
bearing molecules - typical of massi ve hot cores - towards 
IRAS 162293-2422 dCazaux et alJ2003HBottinelli et al.l2004bl) 
and NGC1333-IRAS4A dBottinelli et al.ll2004al) are challeng- 
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ing these theoretical concepts, emphasising that the chemistry 
in low mass protostar envelopes is still poorly understood and 
needs to be further investigated. 

Methanol is particularly suited for the study of hot cores 
around low mass protostars, because it is relatively abundant 
in grain mantles and dominates the organic content of evap- 
oration zones. Indeed, such enrichments of methanol have 
already be en observed in the inn er regions of high mass 
protostars dvan der Tak et alJl2000al). and in the gas shocked 
by p rostellar outflows (Bachiller et al.lll998l buckle & Fullerl 
2002). Moreover, methanol, being a slightly asymmetric-top 
molecule, can be used to probe a large range of physical con- 
ditions, both in density and temperature (e.g. Leurini et al. 
2004). The recent computati ons of collisional r ate coefficients 
of methanol with para-H2 IPottage et alJl2004l) make precise 
determination of these conditions possible. 

In this paper, observations of methanol transitions from 
low mass protostars are presented. A detailed radiative trans- 
fer model is used to derive the methanol abundance inside the 
protostellar envelopes. Finally, the methanol abundances in the 
envelopes are compared to those in other environments, and the 
implications for the formation of this molecule are discussed. 

2. Observations 

Six Class protostars (see Table [0 were observed using the 
Institut de Radio-Astronomie Millimetrique 30 meters tele- 
scope (IRAM-SOm) 1 , and the James Clerk Maxwell Telescope 
(JCMT) 2 . The methanol 5^-4^ transitions were observed with 
the IRAM-30m telescope, while the 7k~6k were observed with 
the JCMT. JCMT and CSO obser vations of IRAS 16293-24 22 
from the literature were also used dvan Dishoeck et al J 19951) . 

The IRAM-30m observations were obtained in November 
2002. The A230 and B230 receivers were used simultaneously 
with the autocorrelator, providing a spectral resolution of 0.5 
km/s and a bandwidth of 320 MHz. The calibration and the 
pointing were regularly checked and were found to be better 
than 20% and 3" respectively. For all the sources the obser- 
vations were done in position switching mode, after checking 
that the reference position was free of line emission. In order 
to determine the spatial extent of the methanol emission, small 
maps of 3 by 3 pixels were made, with a 10" sampling. 

The JCMT observations were obtained in September 2000 
and from February 2001 to February 2003. The B3 receiver 
was used with the digital autocorrelator spectrometer (DAS) in 
setup with a bandwith of 500 MHz and a spectral resolution 
of 0.4 km/s. As for the IRAM-30m observations, the point- 
ing and calibration were regularly checked and was found to 
be better th an 30% and 3" respectivel y. A few lines already 
observed bv Ivan Dishoeck et all i 19951) towards IRAS 16293- 
2422 were re-observed in September 2000. The lines fluxes 

1 IRAM is an international venture supported by INSU/CNRS 
(France), MPG (Germany) and IGN (Spain). 

2 The JCMT is operated by the Joint Astronomy Center in 
Hilo, Hawaii on behalf of the present organizations: The Particle 
Physics and Astronomy Research Council in the United Kingdom, 
the National Research Council of Canada and the Netherlands 
Organization for Scientific Research. 



from the two datasets were compared and were found to be 
in agreement within 30%. The JCMT observations were done 
in beam switching mode with a 180" offset, large enough to 
avoid line contamination by the outflows or by the cloud. 

3. Results 

Fig.^and[2]show the spectra towards the six observed sources. 
Several 5-4 and 7-6 lines are detected in the bands observed 
with IRAM and JCMT respectively. The lines with lowest level 
energy (E up ~ 40 - 70 K) in each band are the brightest in the 
spectra and are detected toward all sources. On the contrary, 
only a few lines with larger upper level energies are detected 
in most sources, with NGC1333-IRAS2 showing the largest 
number of detected methanol lines (26 lines, up to upper level 
energies of 259 K). The relative intensity of the lines with 
low level energy to lines with high level energy varies from 
source to source, suggesting different excitation temperatures. 
The richer methanol spectrum of NGC1333-IRAS2 is, in this 
respect, likely due to a relatively larger excitation temperature 
with respect to other sources. 

The lines with lowest level energy show, in some sources, 
bright high velocity wings, sometimes asymmetric. This 
is the case for NGC1333-IRAS4A, NGC1333-IRAS4B and 
NGC1333-IRAS2. In the other three sources (L1448-MM, 
L1448-N and L1557-MM), all the lines, including those with 
a low level energy, have relatively narrow Gaussian profiles 
(between 1 and 4 kms -1 ) with much weaker broad wings, if 
any. Usually, large asymmetric wings testify to the presence 
of outflowing material. Indeed, the maps obtained at IRAM 
in the direction of NGC1333-IRAS4A and NGC1333-IRAS2 
confirm the presence of large scale outflows probed by the two 
methanol lines with lowest level energy. 

Fig. shows the CH 3 OH 5k-4jc emission map of 
NGC1333-IRAS4A. Towards the central position, the two lines 
with lowest level energy have Gaussian profiles with high ve- 
locity wings. North and south of the source, the wings become 
larger and asymmetric, red to the north-east and blue to the 
south-west. Also, the emission of the two lines with lowest 
level energy is not peaked in the central position, but rather 10" 
south of it. The morphology of the observed emission is con- 
sistent with the direction of the outflow elongating along the 
north-south axis near NGC13 33-IRAS4A^ as seen in CO 3-2 
CS 7 -6 or SiO 2-1 emission JBlake et alJll995t iLefloch et alJ 
ll 99Rh . 

The map of CH 3 OH 5^-4^ emission of NGC1333-IRAS2 
is presented in Fig. |4] In the central position, the lines pro- 
files are similar to those observed in NGC1333-IRAS4A: a nar- 
row Gaussian component, with high velocity wings, although 
less intense. The emission becomes red-shifted north and east- 
ward of the source, and decreases in the west and south di- 
rection. Two outflows have been observed towards NGC1333- 
IRAS2. A large scale outflow, oriented along the north-south 
axis of the source, has been detected in CO 3-2 and 2-1 emis- 

3 There is an abrupt change of the position angle of the CO 3-2 flow 
from approximatively 45° on the large scale to 0° near NGC1333- 
IRAS4A <Blakeetalll99l . 



Table 1. The observed sample. 
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Source 


or(2000) 


(5(2000) 


Cloud 


Dist. a 


I b 


M env b 


L 1 L,\) \ c 


7\,oi c 










(pc) 


(Lq) 


(Mq) 


(%) 


(K) 


NGC1333-IRAS4A 


03:29:10.3 


+31:13:31 


Perseus 


220 


6 


2.3 


5 


34 


NGC1333-IRAS4B 


03:29:12.0 


+31:13:09 


Perseus 


220 


6 


2.0 


3 


36 


NGC1333-IRAS2 


03:28:55.4 


+31:14:35 


Perseus 


220 


16 


1.7 


<, 1 


50 


L1448-MM 


03:25:38.8 


+30:44:05 


Perseus 


220 


5 


0.9 


2 


60 


L1448-N 


03:25:36.3 


+30:45:15 


Perseus 


220 


6 


3.5 


3 


55 


L1157-MM 


20:39:06.2 


+68:02:22 


Isolated 


325 


11 


1.6 


5 


60 


IRAS 16293-2422 


16:32:22.7 


-24:38:32 


p-Ophiuchus 


160 


27 


5.4 


2 


43 



" From Andre et al. 1 2000), except for Perseus sources 1 
b From J0rgensen et al. (2002). 
c From Undre et alJfcOOOh 



sion JKnee & Sandelll EoOO). A s econd outflow, oriented east 
and w est, is detecte d in SiO 2-1 ( Blakjp h 996t [.Tdrgensen et all 
l2004ah and CS 3-2 JLan ger et al]ll996I) . Strong CH 3 OH lines 
have also been observed at the endpoint of this outflow (~ 70" 
east and west of the central position), where the methanol abun- 
dance is enhanced by a factor of 300 (Bac hiller et alll998l) . In 
our map, low energy transitions trace the red lobe (towards the 
north) of the outflow, as well as the blue (towards the west) 
lobe of the east-west outflow, but only narrow and weak lines 
are seen towards the east and south. On the contrary, lines with 
high energy, also seen on the JCMT spectrum, appear only in 
the central position. 

Since this paper focuses on the methanol emission from the 
envelopes surrounding the protostars, in the following we will 
restrain the discussion to this component only and will not an- 
alyze further the outflow component, other than for disentan- 
gling it from the envelope emission. In practice, we separated 
the envelope emission, assumed to give rise to a Gaussian line 
centered on the Vi sr of the source, from the outflow emission, 
assumed to be the residual of the Gaussian. Tables l4l to 1 101 list 
the line fluxes and widths of the Gaussians, for each source, as 
derived from this analysis. As previously said, in most cases the 
lines are rather "clean" Gaussians with narrow widths, and little 
contamination from wings. However, a few lines are strongly 
contaminated by the outflow component, and such a separation 
between the envelope and outflow components was not possi- 
ble. In these cases, a single Gaussian was fitted on the entire 
profile, therefore including both envelope and outflow contri- 
butions. These lines are explicitly marked in Tables I4ltoll0l 

In the three sources where the outflow emission is relatively 
important (NGC1333-IRAS4A, IRAS4B and IRAS2), several 
lines observed with IRAM are significatively narrower than the 
lines observed with the JCMT: between 1 and 4 km s _1 the for- 
mers, and up to 8 km s the latters (before outflow component 
subtraction). Fig.|5]illustrates this difference, and presents two 
methanol line profiles observed with IRAM-30m and JCMT 
towards NGC1333-IRAS2. The line observed with the JCMT 
clearly shows a red high velocity wing, while the IRAM spec- 
trum does not. Likely, this difference is due to the larger beam 
of JCMT (13") with respect to IRAM-30m (10"). The former 
probably picks up more extended emission from the outflows 



than the latter. In addition, or alternatively, this difference can 
be due to a slightly different pointing of the two telescopes. 

The second panel of Fig. [5] shows the case of L1448-N, 
where, on the contrary, the lines from the two telescopes are 
very similar, confirming the lower contribution from the out- 
flow to the line emission in this source. Fig. [5] also show a 
formaldehyde line along with the two methanol IRAM and 
JCMT lines. The figure shows that, in the case of NGC1333- 
IRAS2, the JCMT methanol line is indeed much more contam- 
inated by the outflow than the formaldehyde line, also obtained 
at JCMT. However, inspection of the spectra obtained towards 
L1448-N does not show any substantial difference, supporting 
our conclusion that, in this source, the outflow does not play a 
major role in the methanol emission. 

In summary, profiles and maps of the detected methanol 
lines suggest that in half of the sample sources (NGC1333- 
IRAS4A, NGC1333-IRAS4B and NGC1333-IRAS2) the lines 
with lowest level energy are strongly contaminated by the out- 
flowing gas, whereas in the the lines with the highest level en- 
ergy the problem is much less severe. In the remaining three 
sources (L1448-MM, L1448-N and L1157-MM) no evidence 
of severe outflow contamination has been observed in any 
methanol line. For most of the lines showing contamination by 
the outflow, the outflow emission was disentangled from the en- 
velope emission. For a few lines marked in Tables|4]to[6] such a 
separation was not possible. In all other cases, we are confident 
that the quoted fluxes are representative of the emission from 
the envelope only. 

4. Modeling 

4.1. Rotational diagrams 

Rough estimates of the kinetic temperatures and column den- 
sities have been obtained from the rotational diagrams of each 
source (Fig.|6j. The observations are reasonably well fitted by 
a straight line, with some scattering, which is probably due to 
opacity and non-LTE effects. No differences are observed be- 
tween the JCMT and IRAM-30m observations, despite the dif- 
ferent beam widths of the two instruments. 

The derived rotational temperature and column densities 
are summarized in Table|2] The derived column densities range 
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Fig. 1. On source CH 3 OH 5 k -4-k spectra. 



tional temperatures are around 30 K. Two notable exceptions 
are NGC1333-IRAS2 and IRAS 16293-2422, whose rotational 
temperatures are about 100 and 85 K respectively. That the 
temperatures are larger in these two sources is not a surprise, 
for they are the only two sources where several lines with high 
level energy are detected, despite the fact that lines with low 
level energy are not much brighter than in the other sources. 
It is therefore very likely that the observed methanol lines in 
NGC1333-IRAS2 and IRAS 16293-2422 originate in the hot 
region of the envelope. 



iBuckle & Fuller (2002), methanol can be very sub-thermally 
populated in interstellar conditions, and the rotational temper- 
ature derived from methanol rotational diagram analysis could 
be significantly lower than the kinetic temperature. Part of the 
differences seen from one source to an other could be due to 
differences in excitation. In the next section, a more detailed 
analysis including non-LTE excitation is developed, in order to 
take these effects into account. 



4.2. Jump models 



Rotational temperatures should, however, be considered To determine more precisely the methanol abundance in the 
with some caution, because they may not reflec t the a ctual envelopes, a detailed 1-D radiative transfer model has been de- 
kinetic temperatures. As noted by Bachille r et alJ 1 119951) and veloped. The model uses the escape probability formalism, and 
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Fig. 2. As in Fig.Q]for CH 3 OH 1 K -6 K . 



Table 2. Results from the rotational diagrams. 



Source 



Trot 
(K) 



W(CH 3 0H) 

(cm- 2 ) 



NGC1333-IRAS4A 


24 


± 2 


(5.1 


± 1.0) 


X 


10 14 


NGC1333-IRAS4B 


34 


+ 4 


(3.5 


±0.8) 


X 


10 14 


NGC1333-IRAS2 


101 


± 16 


(3.4 


±0.6) 


X 


10 14 


L1448-MM 


46 ± 10 


(8.8 


±2.4) 


X 


10 13 


L1448-N 


19 


±3 


(9.4 


±3.6) 


X 


10 13 


L1157-MM 


33 


±25 


(5.3 


±4.9) 


X 


10' 3 


IRAS 16293-2422 


84 


+ 8 


(8.1 


±0.9) 


X 


10 14 



has been presented in detail in lPaoer J , The energy levels and 



Einstein coefficients wer e taken from the Co logne Molecular 
Database Spectroscopy JMiiller et all l200l|). CH ^OH colli- 
sional rates with para H2 from IPottag^etal . [ 12004 were used. 
The threefold hindering potential of the molecule leads to the 
formation of A and E symmetry states (see Appendix IA"V 
Because no radiative or collisional transitions are possible be- 
tween the E and A symmetry species, these were modelled sep- 
arately. The first 100 levels were considered for each symme- 
try of the molecule. Finally, only the ground torsional vibration 
state (v f =0) was considered. The first excited state of the tor- 
sional vibration (v,=l) is about 200 cirT 1 (~ 290 K) above the 
ground state. Levels with v,=l might be excited and decay to 
the torsional vibration ground state, but in a higher rotational 
state. However, radiative excitation due to the dust is generally 
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Fig. 3. CH3OH 5 k -4k maps of NGC1333-IRAS4A. The outflow is elongating on a north south axis (see text). 



inefficient, and, for the density and temperature ranges present 
here, c ollisional excitation t o those levels is expected to be neg- 
ligible JSchoier et alJ2002t) . 

The envelopes have been assumed to be spherically sym- 
metric, and power-law density profiles have been adopted. The 
index of the power-law density profiles and the dust tem- 
perature profiles were taken from J0rgensenet al. (2002) and 
IPaoer 1 The methanol abund ance profile has been supposed to 
follow a step function (Paper I). The abundance is set to a con- 
stant value X out , in the outer cold part of the envelope. This 
abundance jumps to a X{ n value in the inner and hotter part of 
the envelope, because of the evaporation of the grain mantles. 

The adopted "jump" abundance profiles may seem simplis- 
tic, because other physical and chemical processes in the en- 
velope could lead to more complex abundance profiles. It has 
been suggested that the formaldehyde abundance in the outer 
envelop e of L1448-C and I RAS 16293-2422 follow a "drop" 
profile ( Scho ier" et al J 12004). In this model, the abundance is 
reduced in the cold dense part of the envelope where CO is 



frozen out, but is undepleted in the inner region of the envelope, 
where T > 50 K. The outer radius of the "drop" region is fixed 
where the density equals 10 5 cm~ 3 . Support for this suggestion 
stems from OVRO observations, which cannot be re produced 
on the short baselines by the jump scenario proposed in Pa per 1 
Increasing the formaldehyde abundance in the outer region was 
found to give a better agreement with the observations. 

Another possible explanation for the discrepancy between 
"jump" models and the OVRO observations could be that, for 
transitions with the lowest level energy, the clouds in which the 
protostars are embedded contribute significantly to the short 
baseline emissio n. T his emission was no t taken into account 
either in PaperJ or in lSchoier et alJ l|2004). Emission from the 
cloud would have the same effect as increasing the formalde- 
hyde abundance in the outer envelope. Indeed, the outer radii 
of the "drop" region is about 5000 AU and 8 000 AU from the 
center for IRAS 16293-2422 and L1448-MM respectively. At 
these radius, the temperature of the gas is only ~ 10 K, and the 
contribution of the cloud is likely to become important. 
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Therefore, the main difference between the "jump" and 
"drop" scenarios is the temperature at wh ich the evaporation 
occurs in the inner region of the envelope. Scho ier et alJ (12004) 
adopted a temp erature o f 50 K, while a temperature of 90 K 
was adopted in IPaper J. Methanol trapped in water ice will 
evaporate at the sublimation temperature of pure water ice 
(~ 90 K). Because observations of interstell ar ices suggests that 
CH3OH i s coexistent with H?Q in the ices ( Ski nner et all 1992: 
iGerakines et ail 1 19991: iBoogert et all EoOO). we have adopted 
90 K for the sublimation temperature. We refer the reader to 
Paper I for a discussion of the effect of varying the evaporation 
temperature. 

In principle, the evaporation temperature depends on the 
composition of ices. If methanol is trapped in water rich (po- 
lar) ices, the evaporation is likely to occur at the evaporation 
temperature of pure water ice (~ 90 K). On the other hand, if 
methanol is trapped in CO rich (apolar) ices, the temperature 
will occur at a lower temperature (~ 30 K). In this study, the 



methanol was assumed to be trapped mostly in water rich ices, 
and the evaporation temper ature wa s therefore assumed to be 
90 K. We refer the reader to Paper I for a discussion of the ef- 
fect of varying the evaporation temperature. 

The values of X out and X m were determined by running a 
grid of 10 x 10 models for each source, with the abundances 
ranging from 10~ 12 to 10~ 4 . The best fit parameters were then 
determined by computing the x 2 at each point of the grid. The 
grid was eventually refined around the best fit solutions to de- 
termine more precisely the abundances. Fig. present the x 1 
maps obtained for each source. The figure indicates the 1, 2 
and 3cr confidence levels, the number of lines considered, as 
well as the reduced^ 2 for the best fit model. The best fit mod- 
els reproduce the observations quite well, with the reduced x 2 
ranging from 0.5 to 3.5. The abundances obtained are presented 
in Table|3] 

The methanol abundance in the outer part of the envelope 
ranges from 3 x 10~ 1() to 2 x 10~ 9 . The inner abundance varies 
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much more from one source to the other. In two sources of 
our sample, NGC1333-IRAS2 and IRAS 16293-2422, the ob- 
servations can only be reproduced by our model if there are 
jumps in the abundances. The inner abundances are 3 x 10" 7 
and 1 x 10~ 7 respectively, i.e. a factor 100 and 200 larger than 
in the outer envelope. We note that IRAS 16293-2422 abun- 
danc es are in good agre ement with those obtained previously 
bv ISchoier etafl ilOOA In NGC1333-IRAS4B and L1448- 
MM, there is weak evidence for abundance jumps (lcr). In 
these sources, the inner abundance jumps to 7 x 10~ 7 and 
5 x 10~ 7 respectively, i.e. a factor about 300 larger than in the 
outer cold part of the envelope. These values should however be 
treated with some caution, because of the low confidence level. 
In NGC1333-IRAS4A, L1448-N and LI 157-MM, the observa- 
tions are well reproduced with a constant CH3OH abundance 
throughout the envelope, even if the presence of a jump up to 
a few 10~ 7 cannot be ruled out by the present observations in 
these sources. 

5. Discussion 

The abundance of methanol in the outer regions of Class pro- 
tostars is 3 - 20 x 10~ 10 , a value comparable with the ones ob- 
served in cold dark clouds (5xl0~ 9 , iDickens et al.l2 000). Also, 
the observations reveal clear evidence for high methanol abun- 
dances in the envelopes of two Class protostars, IRAS 16293- 
2422 and NGC1333-IRAS2. At the two sigma level, two other 
protostars - NGC1333-IRAS4A and LI 448 -MM - also show 
enhanced CH3OH abundances in the inner regions. For the 
other sources, the evidence is inconclusive. 

The presence of large abundance jumps in the inner regions 
of low mass protostars is often cited as evidence for the im- 
portance of evaporation of ice mantles due to heating by the 
newly formed star. These ice mantles are thought to form dur- 
ing the preceding, cold, prestellar core phase. This scenario is 



supported by the detection of high deu terium fr actionati on in 
the warm gas around low mass YSOs jCeccarelli et alJ Fl998: 
lLoinardet al.|l2002l IParise et aljEooi l2004l> . a tell-tale sign 
of low temperature chemistry. The presence of complex or- 
ganic molecules (e.g. CH3OCH3, HCOOCH3) also shows the 
importance of ice evaporation, as these molecules are gener- 
ally thought to form from abund ant ice species through rapid 
reacti ons in the warm gas ( Cha rnlev et alJfl992t ICaselli et alJ 
Il993l) . Methanol plays a central role in this drive towards 
molecular complexity. 

It is clear that the abundance of methanol in the warm in- 
ner regions of low mass protostars can provide important clues 
to the chemistry of regions of star formation and, in particu- 
lar, to the importance of grain surface chemistry versus gas 
phase chemistry. In the remainder of this section, we discuss 
gas phase and grain surface formation routes of methanol and 
compare the abundances derived for the hot core regions with 
models predictions. 

In the gas phase, the main methanol formation route is 
through radiative association of CHJ with H2O, followed by 
dissociative recombination of CH3OH2 (Herbst, private com- 
munication): 



CH+ + H 2 -> CH3OH+ + hv 
CH3OH+ + e~ -> CH3OH + H 



(1) 

(2) 



The bottleneck in this sequence is the first step, the radiative 
association reaction. The rate of this reaction was recently mea- 
sured to be three orders of magnitude low er at 10 K than previ- 
ous estimates. Using the earlier estimate, Robe rts et al.l (E004) 
predicted a methanol abundance in dense interstellar cores of 
10 -8 at early times (10 5 yr), decreasing to 10~ 9 at later times 
(10 7 yr). Decreasing the reaction rate by about three orders of 
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Fig. 6. Rotational diagrams of CH3OH lines. Detected lines are marked by circles with error bars. Arrows indicate 2<x upper 
limits. 



magnitude is expected to decrease the methanol abundance by 
the same factor, i.e. to between 10~ 12 and 10~ n . 

On grain surfaces, methanol is thought to be formed 
by successive h ydrogenation of CO jTielens & HagenlflQsl 
iTielens & Allamandolall987h : 



CO -> HCO -> H 2 CO -> H 2 COH -> CH3OH 



(3) 



This process has been studied experimentally by 
Hida ka et al .] d2004 . and was found to be efficient form- 
ing methanol at low temperature. Theoretical models show that 
methanol is abundant in icy grain mantles (CH^OH/HaO- 0.1) 
when the atomic hydrogen abundanc e in the accreting gas is 
high compared to the CO abundance jKeane & Tielensl 120051) . 

For IRAS 16293-2422 and NGC1333-IRAS2, the methanol 
abundance in the hot core compares well with those found 
in hot cores associated with regions of massive star forma- 
tion (6 - 10 x 10" 8 , ISutton et all Il995t Ivan der Tak et alJ 



2OOOg0) and are much higher than met hanol abundances ob - 
served in cold dark clouds (5 x 10" 9 , iDickens et al.l 12000). 
In contrast, much higher methanol abundances are associated 
with the outflows in the regions of low mass star formation, 
L1157-MM and NGC1333-IRAS2 (2 x 10~ 5 and 2 x 10~ 6 re- 
specti vely, iBachiller & Perez Gutierrez! 1 1997t iBachiller et alJ 
119981) . Likewise, cold ices observed along the line of sight to- 
wards high mass and low mass protostars have methanol abun- 
dances (0.2-2x 10~ 5 ) whi ch are two orders of ma gnitude larger 
than observed in hot cores ( Allamandola et al .ll992HCiibb etaD 
2004; Pontoppida n et alJl2003l l2004 . Lines of sight towards 
background stars on the othe r hand have much lower methanol 
ice abundances (< 5 x 10" 6 , Chia ret al!l996l) . 

It was realized more than a decade ago that the methanol 
abundances observed in the hot core associated with low and 
high mass protostars, interstellar ices and YSO outflows pose 
grave challeng es to gas phase m odels for the formation of 
methanol (e.g. Menten et al. 1988). Adopting the new rate for 
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Table 3. Summary of derived abundances. 



Source 


CO a 


H 2 CO b 


CH 3 OH c 






Outer d 


Inner 6 


Outer d 


Inner 6 


NGC1333-IRAS4A 


8 x 10~ 6 


2 x 10-' u 


2 x 10- s 


7 x 10~ lu 


< 1 x 10~ 8 


NGC1333-IRAS4B 


1 x 10~ 5 


5 x 10-'° 


3 x 10- 6 


2 x 10~ 9 


7 x 10~ 7i 


NGC1333-IRAS2 


2 x 10~ 5 


3 x 10-'° 


2 x 10~ 7 


1 x 10~ 9 


3 x 10" 7 


L1448-MM 


4 x 10- 5 


7 x 10-'° 


6 x 10- 7 


2 x 10~ 9 


5 x 10~ 7i 


L1448-N 




3 x 10-'° 


1 x 10- 6 


7 x 10-'° 


< 4 x 10~ 7 


L1157-MM 


6 x 10~ 6 


8x 10"" 


1 x 10~ 8 


3 x 10- 10 


< 3 x 10~ 8 


IRAS 16293-2422 


3 x 10- 5 


1 x 10- 9 


1 x 10" 7 


1 x 10~ 9 


1 x 10- 7 


L134N d 


1 x 10 4 


2 x 10" 8 


5x 10 


-9 


Orion hot core c 


1 x IQ 4 


7 x IQ- 9 


1 x 10 


-7 


High mass YSO f 


(0.3- 1.7) x 10- 4 


(1 - 10) x 10-' 


(0.4 - 24) x 10~ 9 


(6 - 9) x 10- 8 


LI 157 outflow 8 


1 x 10~ 4 


2 x 10~ 7 


2x 10 


-5 


Ices' 1 


(0.2 - 3) x 10- 5 


(1 - 4) x 10- 6 


(0.2 - 2) x 10- 5 



a From lJOrgensen et alJ l2004bl) and lSchoier et all 1 2002 ). 
b From Pa per 1 
c This study. 

d From lDickens et alJfcOQOh (position C). 
e From lSutton et al J i 1 995TT 

f From lvan der Tak et alJ feOOO aid). 

8 Froml Baehiller & Perez Gutierrez! <1997h . 

h From lEhrenfreund & Charnlev! <2000h and lGibb et all !2004 . assuming X(H 2 0) = 5 x 10~ 5 . 
1 ler confidence level only. 



the radiative association rate, gas phase models fall short by 
four orders of magnitude. Indeed, this is one (additional) rea- 
son why the high methanol abundances in the warm gas are 
generally thought to reflect release into the gas phase of grain 
mantle species in the inner region of the envelope, either due to 
thermal evaporation as a result of the heating by the protostar, 
or due to grain mantle sputtering by the outflow. As discussed 
previously, our maps suggest that the methanol lines with the 
highest level energy do not spatially correspond to the outflow 
emission, but are rather peaked on the central source. Together 
with the small line width, this suggests that this emission origi- 
nates in the hot core where grain mantles evaporate completely 
when the dust temperature exceeds 90 K. A similar conclusion 
was reached f or the hi gh H2CO abundance s observed in these 
same sources (Paper I). Recently, Schoi er et alJ (2004) investi- 
gated the origin of H2CO abundance enhancement, using inter- 
ferometric observations of IRAS 16293-2422 and L1448-MM. 
Their maps reveals a compact emission region and small line 
widths and are not associated with the outflows. They conclude 
that these enhancements originate from th ermal evaporation of 
the grain mantles, as claimed in IPaper 1 Further support for 
the thermal evaporation thesis is provided by the interferomet- 
ric images of complex organic molecules in IRAS 16293-2422, 
which show compact emission around the two protostars of 
the binary system, on s cales smaller than 1" jKuan et all2 004: 
iBottinelli et alJEo04bl . Here, we conclude that the scenario is 
very likely the same for both methanol and formaldehyde. 

If the abundances observed in the hot cores of protostars 
are due to grain mantle evaporation, they should reflect the 
grain mantle composition. The observed discrepancy in the ice 
(0.2-2x1 0~ 5 ) and hot core methanol abundances ( 1 - 7 x 1 0~ 7 ) 
is then somewhat discomforting. Possibly, ice mantles are only 



partially released into the gas phase in these regions or the gas 
phase of hot cores is more heavily processed by gas phase reac- 
tions than realized before JCharnlev et alJl9 92). Such a partial 
release of ice mantles may result from small differences in tem- 
peratures between different size grains (e.g. , smaller grains or 
graphitic grains might be slightly warmer than larger silicate 
grains). Typically, a 10% difference in temperature can cause 
a factor 5 difference in the evaporation rate of ice mantles un- 
der interstellar conditions. Of course, within this scenario, this 
would imply that these warmer grains represent only a small 
fraction of the total ice reservoir. For shocks, partial release of 
the ice mantles may be a natural outcome of the sputtering in 
low velocity shocks. 

The grain mantle scenario can be further tested by con- 
sidering observed abundance ratios, involving CH3OH, H2CO 
and CO rather than absolute abundances. Of course, such a 
comparison inherently assumes that these species evaporate 
simultaneously. For traces of CO, CH3OH and H2CO mixed 
into a F^O-rich ice, evaporation will be domainted by the 
evaporation of the main component. However, some of the 
solid CO may not be coexisting with the H?0-rich ice m antle 
jTielens et alJl99lUChiar et alJl998trPontoppidan et alJ2003l) 
and this rather vola tile ice componen t will evaporate more read- 
ily (Sandford & Allamandola 1993). These uncertainties have 
to be kept in mind in the subsequent discussion. 

The formation of formaldehyde and methanol from CO 
and H accreted ont o grain surfaces has rec ently been stud- 
ied theoretically bv iKeane & Tie lens (2005). In their model, 
the formaldehyde and methanol abundances produced depend 
on two parameters: the density and the ratio P of the prob- 
ability of H to react with CO relative to H2CO (see Eq. [3}. 
Fig. |8] presents the abundances ratio of H2CO and CO rela- 
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tive to CH3OH predicted by the model and observed in var- 
ious environments. First, we note that the H2CO/CH3OH and 
CO/CH3OH abundance ratios observed in the warm gas around 
low mass protostars are somewhat higher than those observed 
in ices (Fig. [HJ. Both sets of observations are well fitted by the 
model, for an adopted P of 1 . The differences between the low 
mass hot cores and the high mass ices reflect a difference in 
the density during accretion, 10 5 versus 10 4 cirT 3 , respectively. 
This value for P is in good agree ment with recent experimen- 
tal estimates faidakaetalJ l2004). We note that for a density 
of 10 5 cm -3 , the calculated CH 3 OH/H 2 ratio is 10 3 . This is 
about a factor of 40 less than the observed CH3OH/H2O ratio 
in IRAS 16293-2422 (a ssuming a water abundance of 3 x 10~ 6 , 
ICeccarelli et aljEoO Oal Higher abundance ratios are obtained 
when accretion occu rs at lower densities wh ere the atomic H 
abundance is higher (Keane & T ielensll2005l) . Indeed, both the 
observed absolute and relative abundances of these species are 
in good agreement with the ice observations towards high mass 



protostar s for a density of 10 4 cm 3 at the time of formation 
ilKeane & Tielen sl2005h . 

We can also compare the observed abundance ratios in low 
mass hot cores with those observed in high mass hot cores 
(Fig. |8j. This comparison reveals that relative to CH3OH, the 
CO abundance is typically enhanced by an order of magni- 
tude while the H2CO abundance is typically less by a factor 
3 in the latter co mpared to the former. As already noted by 
iKeane & Tielensl l|2P05'), the hot core composition in regions 
of massive star formation are quite different from those ob- 
served in the ices towards the same type of objects (cf. Fig. [8). 
Possibly, the composition of the hot core has already been sub- 
stantially altered from that of the evapora ting ices by reactions 
in the warm gas jKeane&Tielensl E005). Hot core chemistry 
models predict that the abundance of H2CO decreases before 
CH3OH ilCharnlevetal]ll992h . so the H 2 CO over CH 3 OH ra- 
tio is expected to decrease with time. Hence, the points in Fig 
|8]are expected to shift from the calculated grain surface chem- 
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towa rds high mass YSO (Ivan der Tak et alJl2000albh . Adapted 
from lKeane & Tielensl J2Q05I) . 



istry curves to the lower right corner with time. From this per- 
spective, the differences in the abundance ratios between high 
and low mass hot cores could therefore be interpreted in terms 
of different ages for these two kind of regions: high mass hot 
cores are typically older than low mass hot core. Gas phase 
chemistry calculations are needed to confirm this scenario. 

Summarizing this discussion, we conclude that current 
gas phase chemistry models fail to reproduce the observed 
methanol abundances by many orders of magnitude. Based 
upon recent experimental and theoretical studies, good agree- 
ment is obtained between the observations and grain surface 
chemistry models in terms of relative abundance ratios of 
CH3OH, H2CO and CO. However, reproduction of absolute 
abundances seems to be more troublesome, perhaps pointing 
towards the presence of multiple ice components and for selec- 
tive partial evaporation. While this comparison presents, at the 
moment, mixed results, further observational studies of these 
and related species, in combination with laboratory studies and 
theoretical models, promises the opportunity of developing a 
chemical clock which can "time" the star formation process. 



6. Conclusions 

Methanol lines observations of a sample of Class protostars 
have been presented. Using a detailed 1-D radiative transfer 
model, the methanol abundance profiles have been derived and 
compared with chemical scenarios of the formation of this 
molecule. Our main conclusions are: 

1. The methanol abundances in the outer cold envelopes (< 
90 K) are 3-20 x 10~ 10 , a comparable value with the ones 
observed in cold dark clouds. 

2. Constant abundance models fail to reproduce the observa- 
tions in two sources of our sample, IRAS 16293-2422 and 
NGC1333-IRAS2. A jump at the radius where the temper- 
ature reaches 90 K has to be included in order to reproduce 
the observations. In two other sources, NGC1333-IRAS4B 
and LI 448 -MM, jumps in the abundance are also detected, 
but at a lower confidence level. In the three remaining 
sources, NGC1333-IRAS4A, L1448-N and L1157-MM, 
the observations are well reproduced by a constant abun- 
dance model, but the presence of a jump in the abundance 
cannot be ruled out. 

3. In the sources where abundances jumps are detected, the 
abundance in the inner hot region (> 90 K) are 1-7 x 10~ 7 , 
i.e. about two orders of magnitude larger than in cold dark 
cloud. The abundances in these hot cores compare well 
with those derived for the Orion hot core, as well as in the 
inner region of high mass YSO. 

4. We have compared the observed methanol abundances with 
gas phase and grain surface chemistry models. Gas phase 
chemistry models predict methanol abundances (10~ n ) 
much lower than observed in these hot cores. Grain sur- 
face chemistry models provide good agreement with the 
observed abundance ratios involving CH3OH, H2CO, and 
CO. However, absolute abundances are more difficult to re- 
produce and may point towards the presence of multiple 
ice components and the effects of difference in volatility 
between different ices species. 
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Table 4. Lines intensities and width for methanol transitions in NGC1333-IRAS4A. 



Transition 


Eup 


V 


J T mb dv 


Av 


Telescope 


HPBW 




(K) 


(GHz) 


(Kkms- 1 ) 


(kms- 1 ) 




(") 


5 -4 E + 


48.0 


241.700 


1.70 ±0.41 


2.7 ± 0.4 


IRAM-30m 


10 


5,-4, E- 


40.5 


241.767 


4.63 ± 0.87 


2.9 ± 0.2 


IRAM-30m 


10 


5o-4 A + 


34.8 


241.791 


7.46 ± 1.40 


3.7 ± 0.2 


IRAM-30m 


10 


5 4 -4 4 A ± 


115.2 


241.806 


<0.21 




IRAM-30m 


10 


5 4 -4 4 E" 


122.8 


241.813 


<0.21 




IRAM-30m 


10 


5 4 -4 4 E + 


135.8 


241.829 


<0.21 




IRAM-30m 


10 


5 3 -4 3 A* 


84.7 


241.832 


<0.21 




IRAM-30m 


10 


5 2 -4 2 A" 3 


72.6 


241.842 


<0.21 




IRAM-30m 


10 


5 3 -4 3 E" 


97.6 


241.852 


<0.21 




IRAM-30m 


10 


5,-4, E+ 


55.9 


241.879 


<0.21 




IRAM-30m 


10 


5 2 -4 2 A + 


72.6 


241.887 


<0.21 




IRAM-30m 


10 


7,-6, E" 


70.6 


338.345 


6.71 ± 2.08 


8.1 ±0.1 d 


JCMT 


13 


7o-6 A +b 


65.1 


338.409 


7.95 ± 2.45 


8.0 ± 0.8 d 


JCMT 


13 


7 6 -6 6 E" 


254.7 


338.431 


<0.06 




JCMT 


13 


7 6 -6 6 A ± 


258.9 


338.442 


<0.06 




JCMT 


13 


7 5 -6 5 E- 


189.2 


338.457 


0.03 ± 0.04 


1.1 ± 1.3 


JCMT 


13 


75-65 E+ 


206.1 


338.475 


<0.06 




JCMT 


13 


75-65 A* 


203.0 


338.486 


<0.06 




JCMT 


13 


7 4 -6 4 E~ 


153.1 


338.504 


<0.06 




JCMT 


13 


7 2 -6 2 A~ c 


102.8 


338.513 


0.46 ± 0.22 


6.9 ± 1.4 d 


JCMT 


13 


7 4 -6 4 E + 


166.0 


338.530 


<0.06 




JCMT 


13 


7 3 -6 3 A* 


114.9 


338.541 


0.75 ± 0.27 


7.4 ± 0.5 d 


JCMT 


13 


7 3 -6 3 E- 


127.9 


338.560 


<0.06 




JCMT 


13 


7 3 -6 3 E + 


112.8 


338.583 


0.14 ± 0.07 


3.1 ±0.9 d 


JCMT 


13 


7,-6, E+ 


86.1 


338.615 


1.40 ± 0.47 


7.6 ± 0.3 d 


JCMT 


13 


7 2 -6 2 A + 


102.8 


338.640 


0.29 ±0.13 


6.5 ± l.l d 


JCMT 


13 



a Blended with 5 3 -4 3 E + . 
b Blended with 7 6 -6 6 E + . 
c Blended with7 4 -6 4 A*. 

d Uncertain flux because of outflow contamination 
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Table 5. As in TableEfor NGC1333-IRAS4B. 



Transition 


Eup 


V 


J T mb dv 


Av 




Telescope 


HPBW 




(K) 


(GHz) 


(Kkms- 1 ) 


(km s 


-') 




(") 


5o-4 E + 


48.0 


241.700 


1.56 + 0.50 


1.2 ± 


0.4 


IRAM-30m 


10 


5,-4, E" 


40.5 


241.767 


2.39 ±0.53 


1.1 ± 


0.1 


IRAM-30m 


10 


5o-4 A + 


34.8 


241.791 


3.74 ± 0.77 


1.3 ± 


0.1 


IRAM-30m 


10 


5 4 -4 4 A* 


115.2 


241.806 


<0.30 






IRAM-30m 


10 


5 4 -4 4 E" 


122.8 


241.813 


<0.30 






IRAM-30m 


10 


5 4 -4 4 E + 


135.8 


241.829 


<0.30 






IRAM-30m 


10 


5 3 -4 3 A* 


84.7 


241.832 


<0.30 






IRAM-30m 


10 


5 2 -4 2 A" 3 


76.8 


241.842 


<0.30 






IRAM-30m 


10 


5 3 -4 3 E" 


97.6 


241.852 


<0.30 






IRAM-30m 


10 


5i-4, E + 


55.9 


241.879 


<0.30 






IRAM-30m 


10 


5 2 -4 2 A + 


72.6 


241.887 


<0.30 






IRAM-30m 


10 


7i-6i E" 


70.6 


338.345 


5.27 ± 1.66 


3.7 ±0.1 


JCMT 


13 


7o-6 A +b 


65.1 


338.409 


5.78 ± 1.81 


3.4 ± 


0.6 


JCMT 


13 


7 6 -6 6 E - 


254.7 


338.431 


<0.08 






JCMT 


13 


7 6 -6 6 A ± 


258.9 


338.442 


<0.08 






JCMT 


13 


7 5 -6 5 E- 


189.2 


338.457 


<0.08 






JCMT 


13 


75-65 E + 


206.1 


338.475 


<0.08 






JCMT 


13 


75-65 A* 


203.0 


338.486 


<0.08 






JCMT 


13 


7 4 -6 4 E" 


153.1 


338.504 


<0.08 






JCMT 


13 


7 2 -6 2 A~ c 


131.2 


338.513 


0.76 ± 0.29 


4.7 ± 


0.5 


JCMT 


13 


7 4 -6 4 E + 


166.0 


338.530 


<0.08 






JCMT 


13 


7 3 -6 3 A ± 


114.9 


338.541 


1.29 ± 0.47 


5.3 ± 


0.4 


JCMT 


13 


7 3 -6 3 E" 


127.9 


338.560 


<0.08 






JCMT 


13 


73-63 E + 


112.8 


338.583 


0.30 ±0.19 


3.0 ± 


1.0 


JCMT 


13 


7i-6, E+ 


86.1 


338.615 


1.60 ±0.61 


4.2 ± 


0.4 


JCMT 


13 


7 2 -6 2 A + 


102.8 


338.640 


0.52 ± 0.27 


4.0 ± 


1.1 


JCMT 


13 



"Blended with 5 3 -4 3 E + . 
b Blended with7 6 -6 6 E + . 
c Blended with 7 4 -6 4 A*. 
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Table 6. As in TableEfor NGC1333-IRAS2. 



Transition 


Eup 


V 


J T mb dv 


Av 


Telescope 


HPBW 






(K) 


(GHz) 


(Kkms- 1 ) 


(kms- 1 ) 




(") 


5n 


-4n E + 


48.0 


241.700 


0.89 + 0.34 


2.9 ± 0.8 


IRAM-30m 


10 


5, 

u 1 


-4, E~ 


40.5 


241.767 


1.67 ± 0.41 


2.4 ± 0.3 


IRAM-30m 


10 


5n 

-'\) 


-4o A + 


34.8 


241.791 


1.74 ± 0.40 


2.0 ± 0.2 


IRAM-30m 


10 


5, 


-4 4 A* 


115.2 


241.806 


0.54 ± 0.20 


3.5 ± 1.0 


IRAM-30m 


10 




-4 4 


122.8 


241.813 


0.47 ± 0.19 


3.3 ± 1.0 


IRAM-30m 


10 


5, 


-4 4 E + 


135.8 


241.829 


0.42 ± 0.22 


3.4 ± 1.3 


IRAM-30m 


10 


5, 

•'J 


-4, A* 


84.7 


241.832 


0.74 ± 0.27 


3.2 ±0.8 


IRAM-30m 


10 


5, 

u 1 


-4o A~ a 

^2 ' 1 


76.8 


241.842 


0.75 ± 0.22 


3.6 ±0.5 


IRAM-30m 


10 


5 3 


-4i E~ 


97.6 


241.852 


0.39 ±0.16 


2.8 ± 0.9 


IRAM-30m 


10 


1 


-4, E + 

-r j j_, 


55.9 


241.879 


1.00 ± 0.26 


3.5 ±0.5 


IRAM-30m 


10 




-Ai A + 


72.6 


241.887 


0.63 ± 0.18 


2.7 ± 0.5 


IRAM-30m 


10 


7, 


-61 E~ 


70.6 


338.345 


1.59 ± 0.59 


5.9 ± 0.5 d 


JCMT 


13 


7n 

' U 


-6n A +b 


65.1 


338.409 


1.89 ± 0.68 


6.5 ± 0.5 d 


JCMT 


13 


7, 


-6 fi E" 


254.7 


338.431 


0.27 ± 0.14 


3.0 ±0.8 


JCMT 


13 


7, 


-6 6 A* 


258.9 


338.442 


0.13 ± 0.12 


3.2 ± 2.4 


JCMT 


13 


7 S 


-65 E" 


189.2 


338.457 


0.27 ±0.19 


3.1 ± 2.0 


JCMT 


13 


7, 


-6<; E + 


206.1 


338.475 


0.27 ± 0.18 


2.7 ± 1.0 


JCMT 


13 


7, 


-65 A* 


203.0 


338.486 


0.29 ± 0.20 


3.0 ± 1.7 


JCMT 


13 


7 4 


-64 E- 


153.1 


338.504 


0.37 ± 0.22 


3.4 ± 1.3 


JCMT 


13 


li 


-6 2 A- c 


131.2 


338.513 


0.40 ± 0.23 


3.2 ± 1.3 


JCMT 


13 


7 4 


-6 4 E + 


166.0 


338.530 


0.43 ± 0.24 


3.2 ± 1.1 


JCMT 


13 


7 3 


-63 A ± 


114.9 


338.541 


0.71 ± 0.34 


4.8 ± 1.0 


JCMT 


13 


7 3 


-63 E~ 


127.9 


338.560 


0.33 ± 0.20 


3.2 ±0.8 


JCMT 


13 


7 3 


-63 E + 


112.8 


338.583 


0.41 ± 0.22 


3.1 ±0.9 


JCMT 


13 


7r 


-61 E + 


86.1 


338.615 


0.67 ± 0.33 


4.2 ± 1.1 


JCMT 


13 


7 2 


-6 2 A + 


102.8 


338.640 


0.37 ±0.19 


2.6 ± 0.7 


JCMT 


13 



"Blended with 5 3 -4 3 E + . 
b Blended with 7 6 -6 6 E + . 
c Blended with 7 4 -64 A*. 

d Uncertain flux because of outflow contamination 



S. Maret et al.: CH 3 OH in low mass protostars, Online Material p 5 
Table 7. As in TableEfor L1448-MM. 



Transition 


Eup 


V 


J T mb dv 


Av 




Telescope 


HPBW 




(K) 


(GHz) 


(Kkms- 1 ) 


(km s 


-') 




(") 


5o-4 E + 


48.0 


241.700 


0.30 + 0.15 


1.3 ± 


0.1 


IRAM-30m 


10 


51-4, E" 


40.5 


241.767 


0.70 + 0.14 


1.2 ± 


0.1 


IRAM-30m 


10 


5o-4 A + 


34.8 


241.791 


0.93 ±0.17 


1.1 ± 


0.1 


IRAM-30m 


10 


5 4 -4 4 A* 


115.2 


241.806 


<0.05 






IRAM-30m 


10 


5 4 -4 4 E" 


122.8 


241.813 


<0.05 






IRAM-30m 


10 


5 4 -4 4 E + 


135.8 


241.829 


<0.05 






IRAM-30m 


10 


5 3 -4 3 A* 


84.7 


241.832 


<0.05 






IRAM-30m 


10 


5 2 -4 2 A" 3 


76.8 


241.842 


0.26 ± 0.07 


3.5 ± 


0.6 


IRAM-30m 


10 


5 3 -4 3 E" 


97.6 


241.852 


0.25 ± 0.07 


2.0 ± 


0.4 


IRAM-30m 


10 


5i-4, E + 


55.9 


241.879 


0.32 ± 0.08 


1.8 ± 


0.3 


IRAM-30m 


10 


5 2 -4 2 A + 


72.6 


241.887 


<0.05 






IRAM-30m 


10 


7i-6i E" 


70.6 


338.345 


0.60 ± 0.36 


7.5 ± 


2.5 


JCMT 


13 


7o-6 A +b 


65.1 


338.409 


0.54 ± 0.24 


5.9 ± 


1.0 


JCMT 


13 


7 6 -6 6 E - 


254.7 


338.431 


< 0.11 






JCMT 


13 


7 6 -6 6 A ± 


258.9 


338.442 


< 0.11 






JCMT 


13 


7 5 -6 5 E- 


189.2 


338.457 


< 0.11 






JCMT 


13 


75-65 E + 


206.1 


338.475 


< 0.11 






JCMT 


13 


75-65 A* 


203.0 


338.486 


< 0.11 






JCMT 


13 


7 4 -6 4 E" 


153.1 


338.504 


< 0.11 






JCMT 


13 


7 2 -6 2 A~ c 


131.2 


338.513 


< 0.11 






JCMT 


13 


7 4 -6 4 E + 


166.0 


338.530 


< 0.11 






JCMT 


13 


7 3 -6 3 A ± 


114.9 


338.541 


< 0.11 






JCMT 


13 


7 3 -6 3 E" 


127.9 


338.560 


< 0.11 






JCMT 


13 


73-63 E + 


112.8 


338.583 


< 0.11 






JCMT 


13 


7i-6, E+ 


86.1 


338.615 


0.30 ±0.19 


7.9 ± 


2.7 


JCMT 


13 


7 2 -6 2 A + 


102.8 


338.640 


< 0.11 






JCMT 


13 



"Blended with 5 3 -4 3 E + . 
b Blended with7 6 -6 6 E + . 
c Blended with 7 4 -6 4 A*. 
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Table 8. As in TableEfor L1448-N. 



Transition 


Eup 


V 


J T mb dv 


Av 


Telescope 


HPBW 




(K) 


(GHz) 


(Kkms- 1 ) 


(km s- 1 ) 




(") 


5o-4 E + 


48.0 


241.700 


0.39 + 0.15 


1.2 ± 0.4 


IRAM-30m 


10 


5,-4, E" 


40.5 


241.767 


1.00 + 0.20 


1.1 ± 0.1 


IRAM-30m 


10 


5o-4 A + 


34.8 


241.791 


1.33 + 0.25 


1.3 ±0.1 


IRAM-30m 


10 


5 4 -4 4 A* 


115.2 


241.806 


<0.05 




IRAM-30m 


10 


5 4 -4 4 E" 


122.8 


241.813 


<0.05 




IRAM-30m 


10 


5 4 -4 4 E + 


135.8 


241.829 


<0.05 




IRAM-30m 


10 


5 3 -4 3 A* 


84.7 


241.832 


<0.05 




IRAM-30m 


10 


5 2 -4 2 A" 3 


76.8 


241.842 


<0.05 




IRAM-30m 


10 


5 3 -4 3 E" 


97.6 


241.852 


<0.05 




IRAM-30m 


10 


5! -4, E+ 


55.9 


241.879 


0.26 ± 0.09 


2.6 ± 0.6 


IRAM-30m 


10 


5 2 -4 2 A + 


72.6 


241.887 


<0.05 




IRAM-30m 


10 


7i-6i E" 


70.6 


338.345 


0.38 ±0.18 


1.8 ±0.3 


JCMT 


13 


7o-6 A +b 


65.1 


338.409 


0.48 ±0.21 


2.0 ± 0.3 


JCMT 


13 


7 6 -6 6 E - 


254.7 


338.431 


<0.06 




JCMT 


13 


7 6 -6 6 A ± 


258.9 


338.442 


<0.06 




JCMT 


13 


7 5 -6 5 E- 


189.2 


338.457 


<0.06 




JCMT 


13 


75-65 E + 


206.1 


338.475 


<0.06 




JCMT 


13 


75-65 A* 


203.0 


338.486 


<0.06 




JCMT 


13 


7 4 -6 4 E" 


153.1 


338.504 


<0.06 




JCMT 


13 


7 2 -6 2 A~ c 


131.2 


338.513 


<0.06 




JCMT 


13 


7 4 -6 4 E + 


166.0 


338.530 


<0.06 




JCMT 


13 


7 3 -6 3 A ± 


114.9 


338.541 


0.11 ±0.10 


2.8 ± 1.6 


JCMT 


13 


7 3 -6 3 E" 


127.9 


338.560 


<0.06 




JCMT 


13 


73-63 E + 


112.8 


338.583 


<0.06 




JCMT 


13 


7i-6, E+ 


86.1 


338.615 


<0.06 




JCMT 


13 


7 2 -6 2 A + 


102.8 


338.640 


<0.06 




JCMT 


13 



"Blended with 5 3 -4 3 E + . 
b Blended with 7 6 -6 6 E + . 
c Blended with 7 4 -6 4 A*. 
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Table 9. As in TableEfor LI 157-MM. 



Transition 


Eup 


V 


J T mb dv 


Av 




Telescope 


HPBW 




(K) 


(GHz) 


(Kkms- 1 ) 


(km s 


-') 




(") 


5o-4 E + 


48.0 


241.700 


<0.07 






IRAM-30m 


10 


5,-4! E" 


40.5 


241.767 


0.51 ±0.13 


1.2 ± 


0.2 


IRAM-30m 


10 


5o-4„ A + 


34.8 


241.791 


0.54 ±0.13 


1.2 ± 


0.2 


IRAM-30m 


10 


5 4 -4 4 A* 


115.2 


241.806 


<0.07 






IRAM-30m 


10 


5 4 -4 4 E- 


122.8 


241.813 


<0.07 






IRAM-30m 


10 


5 4 -4 4 E + 


135.8 


241.829 


<0.07 






IRAM-30m 


10 


5 3 -4 3 A* 


84.7 


241.832 


<0.07 






IRAM-30m 


10 


5 2 -4 2 A" 3 


76.8 


241.842 


<0.07 






IRAM-30m 


10 


5 3 -4 3 E" 


97.6 


241.852 


<0.07 






IRAM-30m 


10 


5,-4, E + 


55.9 


241.879 


<0.07 






IRAM-30m 


10 


5 2 -4 2 A + 


72.6 


241.887 


<0.07 






IRAM-30m 


10 


7,-6, E" 


70.6 


338.345 


<0.21 






JCMT 


13 


7o-6 A +b 


65.1 


338.409 


0.41 ± 0.27 


7.8 ± 


3.7 


JCMT 


13 


7 6 -6 6 E - 


254.7 


338.431 


<0.21 






JCMT 


13 


7 6 -6 6 A ± 


258.9 


338.442 


<0.21 






JCMT 


13 


7 5 -6 5 E- 


189.2 


338.457 


<0.21 






JCMT 


13 


7 5 -6 5 E + 


206.1 


338.475 


<0.21 






JCMT 


13 


75-65 A* 


203.0 


338.486 


<0.21 






JCMT 


13 


7 4 -6 4 E" 


153.1 


338.504 


<0.21 






JCMT 


13 


7 2 -6 2 A~ c 


131.2 


338.513 


<0.21 






JCMT 


13 


7 4 -6 4 E + 


166.0 


338.530 


<0.21 






JCMT 


13 


7 3 -6 3 A ± 


114.9 


338.541 


<0.21 






JCMT 


13 


7 3 -6 3 E" 


127.9 


338.560 


<0.21 






JCMT 


13 


73-63 E + 


112.8 


338.583 


<0.21 






JCMT 


13 


7,-6, E+ 


86.1 


338.615 


<0.21 






JCMT 


13 


7 2 -6 2 A + 


102.8 


338.640 


<0.21 






JCMT 


13 



"Blended with 5 3 -4 3 E + . 
b Blended with7 6 -6 6 E + . 
c Blended with 7 4 -6 4 A*. 



Table 10. As in TableHfor IRAS 16293-2422. 



Transition 


Eup 


V 


J T mb dv 


Av 


Telescope 


HPBW 




(K) 


(GHz) 


(Kkms" 1 ) 


(kms- 1 ) 




(") 


5o-4 E + 


48.0 


241.700 


2.67 ± 0.83 


4.5 ± 1.0 


JCMT 


20 


5,-4, E- 


40.5 


241.767 


4.36 ± 0.97 


3.3 ±0.3 


JCMT 


20 


5o-4 A + 


34.8 


241.791 


4.91 ± 1.03 


3.2 ±0.2 


JCMT 


20 


5 4 -4 4 A* 


115.2 


241.806 


0.72 ± 0.46 


fixed" 


JCMT 


20 


5 4 -4 4 E" 


122.8 


241.813 


<0.10 




JCMT 


20 


5 4 -4 4 E + 


135.8 


241.829 


<0.10 




JCMT 


20 


5 3 -4 3 A* 


84.7 


241.832 


1.99 ± 0.75 


5.3 ± 1.4 


JCMT 


20 


5 2 -4 2 A- b 


72.6 


241.842 


1.97 ± 0.77 


6.3 ± 1.7 


JCMT 


20 


5 3 -4 3 E" 


97.6 


241.852 


0.61 ± 0.44 


fixed" 


JCMT 


20 


5,-4, E + 


55.9 


241.879 


2.29 ± 0.78 


4.9 ± 1.2 


JCMT 


20 


5 2 -4 2 A + 


72.6 


241.887 


1.67 ± 0.82 


7.8 ± 3.0 


JCMT 


20 


5 2 -4 2 E" 


60.8 


241.904 


3.28 ±0.91 


4.6 ± 0.7 


JCMT 


20 



" Gaussian fit with a fixed line width of 5.0 km.s 
b Blended with 5 3 -4 3 E + . 
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Appendix A: Spectral designation 

Methanol (CH 3 OH) is a slightly asymmetric rotor with hindered internal rotation of the methyl (CH 3 ) group. The first excited state of the 
torsio nal vibration is about 200 cirT 1 above the ground state, low enough to be populated in some particular interstellar conditions I Men ten et alJ 
1986). Energy levels are labelled by the total angular momentum J, its projection along the symmetry a-axis k or K = \k\, the torsional vibration 
quantum number v, and the torsional sublevel quantum number a. The threefold symmetric torsional barrier in methanol inde ed leads to the 
existe nce of three torsional symmetry states, usually designated A, Ei and £2 corresponding to cr = 0, +1,-1, respectively iLees & Baked 
1968). For the E species, k > states are degenerate with the £2 k < st ates, and vic e versa. It is thus simplest to refer to these states as doubly 
degenerate states of E symmetry with positive and negative values for k <Leesll973l) . Levels of the A species are torsionally doubly degenerate 
but this degeneracy is lifted by the slight asymmetry of the molecule. These doublets may be labelled by the C3,. sy mmetry group notatio n Ai 
and A 2 but the most common convention is tu use the labels A + and A~. Full details about notations can be found in lLees & B aker 1 1968J). An 
alternate convention for labelling the A symmetry states is the use of the standard asymmetric top notation Jk„k c , where K a denotes the above 
K, as done for example in the Cologne database for molecular spectroscopy iMiill eret all200ll) . In this paper, the former notation is employed. 



